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Abstract 

We describe an experimental search for spontaneous formation of flux lines 
during a rapid quench of thin YBa2Cu307_5 films through Tc- This effect is 
expected according to the Kibble-Zurek mechanism of a creation of topological 
defects of the order parameter during a symmetry breaking phase transition. 
Spontaneously formed vortices were previously observed in superfluid ^He, 
while a similar experiment in superfluid '^He gave negative results. Using a 
high Tc SQUID, we measured both the magnetic flux in the sample during a 
quench with a sensitivity of 20 (po/cia^, and the field noise which one would 
expect from flux lines pinned in the film. The sensitivity was sufficient to 
detect spontaneous flux at a level corresponding to 10^^ of the prediction. 
Within our resolution, we saw no evidence for this effect. 

PACS:74.40.+k, 67.40.Vs 
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I. INTRODUCTION 



If a system undergoing a phase transition into an ordered state is quenched through 
the phase transition fast enough, topological defects can be created due to the evolution 
of uncorrelated regions of the newly formed phase, having different values of the order 
parameter. The defects appear at the boundary separating several coalesced regions of 
this kind. Such a scenario was first proposed by Kibble in the context of the Grand 
Unified Theory, in order to describe the symmetry-breaking phase transition in the early 
universe. Zurek developed this idea to predict the initial density of defects created during 
the phase transition and suggested specific experiments on condensed matter systems to 
test this scenario. Because of the generality of the theory of phase transitions, topological 
defect formation should occur in every physical system having a relevant symmetry breaking 
of the order parameter during the transition. The main arguments used by Zurek rely 
on the critical slowing down of fluctuations and the divergence of the coherence length 
near a second order phase-transition. Both quantities are influenced dynamically by the 
characteristic time required to complete the transition. In the last few years, experiments 
were carried out on several systems: first, nematic liquid crystals undergoing an isotropic- 
nematic transition |^,^ (In this case the topological defects are disclinations). Another 
system is liquid ''He crossing the A transition as a result of rapid depressurization [^0; 
in this case the topological defects are quantized vortex lines. A third system is liquid 
^He undergoing its superfiuid phase transition In this experiment, the quantized 

vortices are formed during a thermal quench induced following an exothermic neutron- 
induced nuclear reaction. The experiments on liquid crystals and '^He gave results consistent 
with Zurek's prediction. In contrast, after early claims of observing this effect in liquid ^He 
0, a recent improved experiment showed no spontaneous nucleation of vortices within the 
experimental resolution 0. Here, we report an analogous experiment with a superconductor; 
the topological defects are quantized flux lines. Specifically, the experiment aims to observe 
spontaneous flux lines generated during thermal quench of YBa2Cu307_5 (YBCO) thin 
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films through the normal to superconductor transition. The additional importance of this 
experiment is to test the " cosmological" scenario of Kibble and Zurek [Q,^ in a system with 
a local gauge symmetry, where the theory is less clear then in a systems having a global 



gauge (e.g. He and He), due to the evolving gauge field (B) during the transition ||TO|,|TT 



Within our resolution, we found no spontaneous vortex formation down to a level of 10 ^ 
of the predictions, which is in variance with the original theory of Zurek. 



II. THE KIBBLE-ZUREK MECHANISM IN A SUPERCONDUCTOR 

Flux lines may become created spontaneously in type II superconductors during a rapid 
quench through Tc. The most practical way is a thermal quench (which is more reliable 
from pressure quench in these materials, since the pressure dependence of is rather weak). 
Low temperature superconductor (LTS) have a second order phase transition which is well 
described by the Landau-Ginzburg theory and a rather small critical region, thus, the an- 
ticipated initial flux line density should be well predicted by Zurek's theory 0. According 
to these predictions, the initial vortex density after the quench should be: 

^ Uro/rgY (1) 

^0 

Here, and tq are the coherence length and the relaxation time of the order parameter 
at T = respectively. The typical quench time, tq, is defined as: tq = (^^/j^) _^ , £ being the 
reduced temperature. The exponent p is related to the critical exponents of the coherence 
length and the relaxation time as: p = {-^), where: ^ = l^l''^ , t = Tq kP^- From an 
experimental point of view, the problem with LTS is that their coherence lengths are of 
the order of lOOnm and therefore the fiux line density predicted by Eq. (||) is quite small. 
A more favorable situation exist in high temperature superconductors (HTS) which have 
a much shorter coherence length (^o — 2nm). The phase transition of HTS is closer to 
the 3D-XY model, and has a wider critical region than LTS, but as both the coherence 
length and the relaxation time diverge at that should be sufficient for spontaneous 
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vortex formation (provided the quench is fast enough). Another important aspect unique 
to a superconductor is flux pinning, which can cause a significant reduction of the mutual 
annihilation of vortices and antivortices generated during the quench. Pinning can also 
prevent the flux from being expelled out of the film. Finally, within the last few years it 
was found that the order parameter in HTS has a predominant d-wave symmetry, while in 
metal-LTS it is always s-wave. This type of paring can lead to spontaneous flux generation in 
some special configurations (ref. [0 and other ref. therein), but in a homogenous material it 
should not affect the Kibble-Zurek scenario. The values pertinent to YBCO are — 1-6 nm 
, To ^ 5 ■ 10^^^ sec [Q. In our experimental setup, described below, we achieve a 20K/sec 
cooling rate, giving tq ~ 5 sec . 

Because the phase transition in YBCO is not a pure Ginzburg-Landau type, closer to 
the 3D-XY model, the value of the exponent p is deduced from experiments that measured 
the exponents u and /i ||T5HT8|. We chose average values: u ~ 0.67, /i ~ 3.4 (note that this 



is an unusual scaling), and then p ~ 0.3. Due to the short coherence length, the predicted 
initial flux line density generated in the film by a thermal quench is very large: rii ~ 10^° 
cm~^ (vortices and antivortices). An important quantity is the reduced temperature Tat 
which the system returns to the usual critical behavior during the quench, after an initial 
period during which it is out of equilibrium. Spontaneous flux nucleation takes place around 
this temperature, and thus determines the initial vortex density. According to the theory: 
^ — (tq/tq)^^*'^^^^- In our experiment e' is of the same order of magnitude as the width 
of the regime in which strong fluctuations can erase flux lines {sf ^ 10^^), as deduced 
from measurements of flux line noise in thermal equilibrium [^. However, as soon as 
the temperature decreases and the flux lines become pinned, one can hope to preserve a 
significant fraction of the initial flux line density. In our system, the situation is more 
favorable in this respect than in superfluid ^He, since there the regime of strong fluctuations 
is wider. Although the value of e is not necessarily determined by Eq from the Ginzburg 
criterion |^ (in ^He: Eg ~ 0.2) it is not reasonable that it can be smaller by more than 
a few orders of magnitude. In the recent experiment ? ~ 10~^ . Moreover, there is no 



pinning of vortices in bulk ^He which can help to preserve a significant fraction of these 
vortices. The freeze-out coherence length ^, corresponding to is ^ = ^o{t~q/tqY/'^^^^^\ We 
get ^ ~ O.lyum, of the same order of magnitude as the thickness of the films (50-300nm). 
Therefore, the initial vortex array is two dimensional (but the physical system is 3D). In 
our experiment we can measure directly the difference between the number of vortices and 
antivortices, namely the net flux. If the simple picture of well deflned phase regions with 
separation of ^ and with a choice of a minimal phase gradient between these regions (the 
geodesic rule) being strictly correct, than the r.m.s net-flux should be ~ n'J^ . However, 
these arguments are unlikely to be strictly correct, especially in a real superconductor, due 
to the local gauge symmetry with the evolving magnetic fleld ||T0Hl2[. Moreover, the net- 
flux cannot be determined from the simulations done until now p2|,^ because these were 
performed using the simple geodesic rule and flxed periodic boundary conditions. For a 
plausible estimate of the net-flux one must allow some relaxation of the geodesic rule [jl2| . 



This idea was proposed for flrst-order bubble collision, but it can be applied in principle 
to Zurek's theory at the stage at which uncorrelated phase regions with separation of ^ are 
considered. If, for example, we relax the geodesic rule in such a way as to allow a variance 
of one random flux line in an area of 100 (which may contains 100 vortices, antivortices 
and homogenous sites), the r.m.s. net-flux will be determined through a random walk count 
of rij/lOO flux lines. In this case, the net flux is ~ ]^^A'^ .When relaxing the geodesic rule, 
we can count the net topological charge by integrating around the boundaries of the system, 
only after we identify the vortices inside, as otherwise it will not be self-consistent. Under 
this assumption, one sums the phase differences around any given loop while allowing, with 
some probability, a gradient of the phase between adjacent regions which is not minimal. 
The total topological charge will depend on this probability, and may acquire different values 
for the same initial arrangement of the regions having different values of the order parameter. 
One can see that even when almost does not change (a small relaxation of the geodesic 
rule), the value of the net-flux can become signiflcant and rapidly approach it's maximal 
statistical value of ^JrLl. Based on the above assumptions, we estimate the realistic order of 



magnitude of the net flux line density as ^^y^ = lO^cm"^. Obviously, the net density of flux 
is not affected by mutual annihilation, since an equal numbers of vortices and antivortices 
disappear. It may decrease, however, as a result of flux lines being expelled out of the film. 
We show below that under the conditions of our experiment most of the net-flux should 
remain in the film. Our setup, described below, can resolve net flux down to a limit of 20 
(po/cw?. Note that this sensitivity, with the assumption about the net flux being ^y^, 
is sufficient to see the effect even if rii is 4 orders of magnitude smaller than predicted by 
Eq.(|l]). The above estimation was done for a homogeneous phase transition. In our samples 
the superconducting transition is not completely homogeneous for two reasons: temperature 
gradients arising during the quench and a slightly different in different regions of the film. 
We estimate the maximum temperature gradients as VT ~lK/cm. The spread of in the 
different regions of the film is of the same order. The homogeneous approximation holds if 
vt > 's P,p^, where vt is the velocity of the phase transition front propagating across the 



film and 'sis the characteristic speed at which superconducting order parameter fluctuations 
propagate at e". This condition is identical to imposing the demand that VT < TcE/^. The 
estimated value of the RHS of the inequality is more then 10'^ K/cm, and we therefore 
conclude that the homogeneous approximation is correct. 

One more possible experimental geometry is a superconducting loop. One can try to 
measure spontaneous flux generated during a rapid thermal quench of this loop. The basic 
idea here is the same, namely that uncorrelated regions will be generated during the quench, 
with random phases of the order parameter in each one of them. The accumulated phase 
difference around the loop will create a supercurrent and magnetic flux through the center of 
the loop. If the width d of superconductor forming the loop is of the order of magnitude of 



^, the average (rms) number of flux quanta generated in the loop is: ra^ = 1/4-^/ L/^, where 
L is the circumference of the loop 0. In our typical pattern, L ~ 20mm, and d ^ 10/xm 
(somewhat bigger then C.) so we can substitute d instead of ^ and get ~ 11. This simple 
argument is subject to the same uncertainties as the ones encountered above for bulk film, 
due to the presence of the magnetic field. 



III. EXPERIMENTAL SETUP 



In order to measure the flux, we used a high SQUID placed close to the superconduct- 
ing fllm. The SQUID can detect (a) the net-flux nucleated as the film is quenched through 
Tc, and (b) the field noise caused by random hopping of these fiux lines. The sign of the 
net flux should be random for each individual quench. In contrast, the r.m.s. power density 
of the noise spectrum is the same for vortices and antivortices, and should characterize the 
total density of flux, rather than the net density. 

The essential part of the experimental setup is shown in Fig. 1. It is basically a cryostat 
divided into two cells separated by a thin mylar barrier. The SQUID and the sample are 
mounted on the two sides of the barrier, facing one another. The distance between the 
SQUID and the superconducting fllm is about 1mm. The cryostat is immersed in liquid 
nitrogen at 77K. The whole system is carefully shielded from the earth's magnetic held by 
several /x-metal layers arranged inside a soft iron container. The residual magnetic fleld in 
the cryostat was less then 0.2mG. Additional small coil adjacent to the sample was used 
to null this fleld, as well as for testing the fleld dependence of the results. The reason for 
using two separate chambers is that in this arrangement the SQUID remains at a constant 
temperature, while the sample can be heated and cooled independently. To avoid spurious 
magnetic flelds generated by the current used in resistive heating, the fllm is heated using a 
focused light beam. The light is introduced into the cryostat via a quartz rod, terminating 
at about 2mm from the sample. The light illuminates the whole sample area (lOxlOmm) 
and is conflned within the plastic holder tube of the sample. In order to achieve maximum 
efficiency of the heating, the superconducting film is coated by a graphite layer which absorbs 
the fight. 

Cooling of the sample after the heating stops is via a strong thermal link to liquid 
nitrogen, through helium exchange gas present in the cell. The cooling rate through Tc can 
be regulated by changing the pressure of the gas, and hence its thermal conductivity. It 
takes about 1-2 sec to heat the sample above Tc (~90K). Cooling begins immediately when 
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the light is turned off (by a shutter at the top of the cryostat) . The maximum coohng rate 
through T(, is 20K/sec . Because the heating and coohng is done mainly in perpendicular to 
the plane of the substrate, the temperature of the sample is approximately the same along 
its lateral dimensions. Moreover, the critical slowing down of the fluctuation in the film 
near Tc does not affect the coohng rate. It is possible to keep the SQUID's side of the cell 
in vacuum to avoid any heat leak from the sample to the SQUID, in practice, it was not 
necessary. 

To measure the temperature of the sample in real time during the heating-cooling cycle 
we used a thin graphite strip painted on the back side of the substrate. It's heat capac- 
ity is small, similar to that of the superconducting film, and the changes of its resistance 
with temperature enable us to measure the true temperature of the film in real time. The 
temperature dependence of the resistance of the graphite was calibrated against a diode 
thermometer. A typical temperature measurement during the quench is shown in Fig. 2. 
During the heating period, the graphite strip may be hotter then the film by about IK, 
because it faces the light directly, but the cooling rate is approximately the same. We usu- 
ally raised the temperature of the graphite to about lOOK to ensure that all the film was 
heated above T^. We ascertained that the small measuring current through the graphite strip 
did not affect the output of the SQUID, which was the same with or without temperature 
measurement. 

At different times through the experiment, we used two kinds of HTS DC-SQUIDs: a 
commercial M-2700 unit made by Conductus, and a grain boundary SQUID made in our lab. 
Both were operated in a flux locked loop (FLL). The commercial SQUID has a magnetic field 
sensitivity of 2 • 10~^ Gauss/0o and an integrated pick up loop with an area close to Icm^. 
The home made SQUID, made from a thin YBCO film on a SrTiOs bi-crystal substrate, has 
a similar field sensitivity, however, its noise is larger by a factor of 30 than the commercial 
SQUID. The advantage of using the home made SQUID is that it can be mounted closer 
to the sample, as it is not encapsulated. The distance between the home made SQUID and 
the YBCO film is about 0.5mm, while for the commercial unit it is 1mm. Except for the 
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aforementioned S/N ratio, the results presented below were the same with both SQUIDs. 

Working with a HTS SQUID enabled us to keep the whole experiment at 77K and to 
perform the fast thermal quench as described. Measurement of net flux was done contin- 
uously during the heating and cooling cycle, but the noise becomes well defined just close 
to equilibrium (i.e., at the end of the cooling). If a net flux is generated in the film during 
the quench it should be seen as an offset relative to the zero flux state prevailing when the 
film is above Tc. As stated above, according to the Zurek scenario this flux should increase 
in amplitude with the cooling rate while it's sign should be random from one quench to the 
next. If a large number of vortices and antivortices is preserved in the film at the end of the 
cooldown to 77K, the residual flux noise should be stronger compared to a situation when 
the film is relatively clean from vortices (for example, after a very slow cooling or when the 
temperature is above Tc). 

We now turn to estimate the sensitivity of the experiment. The coupling ratio between 
the pick up loop and the SQUID is 10^ and the DC-fiux level change which one can clearly 
resolve is ~ (po/lOO. According to our measurements of the screening of the SQUID by the 
sample, a single vortex in the sample induces (on average) about O.50o into the pick up 
loop. Dividing by the coupling ratio we get a net flux resolution (in the sample) of ~ 20 
net-0o/cm^. This measured value of the coupling coefficient is consistent with calculations 
19| , |20[| . In Fig. 3. we show an example of a sequence in which an external field was applied 



while the sample was cold. The field change picked up by the SQUID is about half of that 
applied, which serves to estimate the coupling coefficient. As the film is heated above T^, 
the whole field penetrates the sample and is picked up by the SQUID. We found the same 
coupling coefficient by following the reverse procedure, namely cooling the film in a field, 
then turning off the external magnetic field, and finally releasing the remanent flux pinned 
in the sample by heating it above T^. From these we conclude that the samples have strong 
flux pinning and that the coupling coefficient of flux into the SQUID is about 0.5. We used 
the same method to test the superconducting continuity of the ring patterns. In addition, 
we found that external magnetic fields of up to lOmG become trapped during cooldown in 



the film very close to Tc- Thus, any spontaneously generated flux remaining in the fllm after 
it is cooled to a temperature outside the interval e below should survive the cooldown 
to 77K. A net flux line density of lO^0o/cm^ is equivalent to a fleld of 2mG. The pinning 
site density in similar fllms was estimated in ref. |^0[ (and in ref. therein) as 1-6-10^° cm~^. 
Thus, we conclude that most of the net flux generated during a quench should remain pinned 
in the fllm during the cooldown. 

To estimate the expected fleld noise in the fllm, we rely on noise measurements in YBCO 
thin fllms, done with a SQUID having fleld sensitivity similar to ours It was found 

that the noise power spectrum is linearly proportional to the magnetic fleld in the fllm and 
has a 1// dependence. The magnitude of the r.m.s. fleld noise (at lOHz) with 10 Gauss 
(lO^vortices/cm^) was 10^^ — 10^^ Gauss/ V Hz at 77K (the exact number depending on the 
sample). The fleld noise of our SQUID is ^ 3 ■ 10^^ Gauss/ V Hz. Thus, we expect to resolve 
a contribution to the noise in the 1// regime for a total flux line density equal or larger 
than 10^ cm~^. This would be possible if 10~^ (or more) of the initial density survives. The 
estimates done so far refer to unpatterned fllms (Icmxlcm area). The limit of sensitivity 
in experiments on fllms patterned into loops (with diameter of ~ 7mm) is ~ 10 0o net. 

Finally, the residual magnetic fleld in the experimental cell (before nulling), was found 
to be less then 0.2mG. Since any nulling procedure is imperfect, a small constant fleld is 
always present in the fllm when cooling through T^.. In order to evaluate the effect of any 
residual external fleld on our data, quenches were performed under different flelds in the mG 
range, varying systematically both in magnitude and in the sign. 



IV. RESULTS AND DISCUSSION 



In this study, we used several types of c-axis oriented epitaxial YBCO fllms: (1) DC- 
sputtered YBCO on (100) SrTiOg; (2) DC-sputtered YBCO on (lOO)MgO; (3) Laser ablation 
deposited- YBCO on (100)SrTiO3; and (4) DC-sputtered YBCO on (001)NdGaO3. We 
have also tested a sintered YBCO ceramic sample. There are several important differences 
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between the various types of films. Generally, sputtered films have better crystallinity than 
the ones grown by laser ablation. The degree of crystallinity is a measure of the defect density 
in the film, and thus the density of the pinning sites. Films grown on (001)NdGaO3 are 
unidirectionally twinned, while those grown on other substrates are bidirectionally twinned. 
In addition, the films on (001)NdGaO3 are lined by unidirectional nano-cracks, perpendicular 
to the twins, and can be thought of as composed of superconducting strips several fim wide. 



separated by the nano-cracks p5[. Consequently, in those films all the flux lines are very 
close to some boundary and can be therefore expelled from the film during cooldown. Films 
of each type, and of varying thickness in the range 50nm-300nm, were tested both as grown 
( 10x10 mm^) and after patterning. The patterns used were: (1) discs (8mm diameter); (2) 
rings of 7mm diameter and 10 to 100/xm width; (3) an array of 100 rings of 200/im diameter 
and 3/im width. 

An experiment is taken by performing a thermal quench as described while recording the 
SQUID's-FLL voltage output vs. time, or equivalently, vs. temperature. For each sample 
we recorded many consecutive quenches. Typical examples of such measurements performed 
on three different kinds of samples are shown in Fig. 4, along with a reference measurement, 
that of a substrate coated with graphite, but without any superconducting film. We show 
here two consecutive quenches for each sample. It is clearly seen that there are small flux 
jumps during the heating and cooling, but no residual flux after the samples are cooled. 
These flux jumps may be different from sample to sample. However, for any given sample 
the flux jumps seen here appear the same through all the quenches, both in the sign and 
magnitude. Moreover, these flux jumps are not affected by external magnetic fields or from 
the field present at the feedback-modulation coil of the SQUID. Thus, these jumps are 
definitely not related to the Kibble-Zurek mechanism. Possible origin of this effect, which 
occurs near T^., will be discussed elsewhere. The fact that we see no residual flux at the end 
of the cooling shows that within our resolution no flux lines were spontaneously nucleated in 
the film during the quench. Our estimate of the predicted net spontaneous flux line density 
is ~ 10^ cm~^. This should give a magnetic field of ~ 2 mG at the sample (ImG at the 
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SQUID, which is about 5-10^ times bigger then our resolution). If spontaneous flux of this 
magnitude was created, we should have seen a large signal during the cooling (about a factor 
of 100 larger than the reproducible flux jumps shown in Fig. 4). Even if the annihilation 
of vortices and antivortices and the expelling of flux out of the film is very fast, we should 
have still seen a transient signal of varying sign and magnitude as the film cooled through 
Tc. Such signals were not observed even at the fastest recording rate ( lO^data points/sec, 
corresponding to a temperature change of 2xl0~^ in e between successive points) in the 
vicinity of Tc. At this sampling rate, the time tm between successive data points is two 
orders of magnitude smaller than the time t = (toTq)^^*^^^^-* during which the system cools 
from Tc to ? during the quench. In our system, t ~ lO^^sec (the recording rate was limited 
by the minimal integration time constant of the Flux Locked Loop). Similarly, negative 
results were found also with the films patterned into rings. Finally, we did not observe any 
increase in the magnitude of the noise at the end of the cooldown following a quench. 

In conclusion, in our experiments we found no evidence for spontaneous flux line forma- 
tion down to a limit of, at least, 10^^ of the predicted initial density. One possible reason 
for the negative result is that the relation for (Eq. (|T])) is just a first estimation and in 
practice it may be smaller. For example, results obtained in experiments on ^He and ^He 
and also in numerical simulations are smaller by as much as one or two orders of magnitude 
than predicted by Eq. (|I]). The negative results in ^He may be a consequence of a fast 
decay of vortex loops which takes place before the measurement even begins. This scenario 
follows from the calculations done by Williams [^]. As he pointed out, a similar mechanism 
may also cause a decay of vortex loops in a superconductor. In our case, vortex loop can 
be thought of as a vortex-antivortex pair. Since this scenario involves only loops, no initial 
imbalance between the number of vortices and antivortices is permitted. If however, vortices 
and antivortices are not created in equal numbers (as loops), then the imbalance between the 
number of vortices and antivortices (the net flux) cannot decay via this mechanism. There is 
no restriction within the Kibble-Zurek picture that the number of vortices and antivortices 
should be equal. In our experiment, the mechanism proposed by Williams would be reflected 
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in the decrease of the total vortex density and hence in the residual noise level at the end 
of the quench, but it will not affect the anticipated net flux density (our main measuring 
technique). The theory of Zurek implies that after the initial creation, the density of vortices 
will decay as ~ riit/t , where t = at Tc. In the ^He experiment, the measurement began 
at a time which is much bigger then t, while in our fast measurements we took data with a 
time resolution t„i <^ t, and we covered the whole interval from t < up to t ^ t. Thus, 
we should have seen something at least during the initial stage, disrespectful of the details 
of any subsequent decay. In addition to the work of Williams, a reexamination of the Zurek 



scenario for the '^He case was done by Karra and Rivers Their calculation yields a 

vortex density consistent with the negative experimental result. It is not clear what should 
be the result of a similar calculation for a 3D high superconductor. A very interesting 
simulation of a rapid cooling of a normal spot created inside a 2D superconductor (including 
pinning) was carried out recently [|28|. In particular, this simulation shows the great im- 
portance of pinning in preserving a significant fraction of the vortex population originating 
in the quench, essentially ad infinitum. It would be highly desirable to extend this work 
to the case where the topological charge during the quench is not conserved. Finally, we 
point out that the very question of two separate superconductors having a well defined phase 
difference before coming into contact is still under debate We hope that our results 

may shed some light on this interesting issue. 
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Figure captions 

Fig. 1. A schematic layout of the essential part of the experimental system. For clarity, 
only the edge of the mylar barrier is shown. The sample itself is held inside a light-tight 
plastic holder to prevent the light from reaching the SQUID. The assembly shown is 
surrounded by several mu-metal shields. 

Fig. 2. A typical profile of the temperature vs. time during a heating-coohng cycle. The 

inset shows an expanded view of the cooling profile near Tc. 

Fig. 3. Response of the system to an application of an external field of 0.3 mG. Initially, 
the film is at 77K, and after the heating-cooling cycle it cools back to this temperature. 

Fig. 4. Typical time dependence of the flux picked up by the SQUID for several samples 
undergoing heating-coohng cycles. The vertical line in the middle of the figure sepa- 
rates two consecutive heating-cooling cycles. The heating interval takes 1-2 sec at the 
beginning of each cycle. Note that the cooling rates for each sample are different, as 
shown by the different times on the respective horizontal axes. 
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